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The S2 site of the KcsA K"*" channel has eight carbonyl ligands in the ion-binding site. A recent 
study suggests that the K"''-over-Na+ selectivity of the S2 site can be understood by noting the larger 
free energy change involved in enforcing an eight water coordination state around Na+ relative to 
K+. The free energies were obtained from the probabilities of observing eight water molecules 
within a coordination sphere whose radius (Ai^jj+ ~ 3.1 A and Aj^+ ~ 3.5 A) extends to the first 
minimum of the ion-water oxygen pair correlation function. Curiously, using the same coordination 
radius led to results that question the very idea of using coordination states in water to understand 
selectivity in the channel. We show that density fluctuations in neat water at the length scale of the 
coordination volume enter the description of the hydration thermodynamics of the ion expressed in 
terms of its coordination states. Density fluctuations explain the sensitivity to the choice of radius. 
After accounting for this effect, the results lead to the conclusion that free energy changes involved 
in the transition between coordination states in water are inadequate to explain selectivity in the 
ion-channel. 



I. INTRODUCTION 

Molecular simulations of the selectivity of the S2 site 
in the KcsA K+ ion channel show that [H, 0, Q 



= AAi'=''(S2) - A/i°^(aq) 
K, 6 kcal/mol. 



[AiNa+(aq)-AiK+(aq)] 



(1) 



Here /S.fjp'^ is the free energy of transferring Na+ from the 
bulk to the S2 site relative to the same quantity for K+. 
/i™ is the excess chemical potential of the ion; /j,^(aq) 
(X = K'^jNa'^) is the hydration free energy and A'^(S) 
is the analogous quantity in the S2 site. 



A. Inferring selectivity from the distribution of 
coordination states of the aqua-ions 

In an effort to rationalize the selectivity of the chan- 
nel, Bostick and Brooks [1] suggest that the selectivity of 
the S2 site can be understood based on the probabilities 
of observing specified n-coordinate structures (X[H20]n) 
of the aqua-ion in bulk water. (The identification of a 
particular n-cluster requires the specification of a coor- 
dination sphere around the aqua-ion.) The idea is that 
the lower probability of observing the eight-coordinate 
state of Na+(aq) relative to K+(aq) indicates a higher 
free energy change involved in enforcing an eight coordi- 
nate state around Na"'"(aq) relative to K+(aq). But the 
S2 site also provides an eight-ligand site, albeit the lig- 
and is a carbonyl and not water. But the behavior in 
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water itself suggests that the "control of the permeant 
ion's coordination state" H is the basis of selectivity. 

Bostick and Brooks [4| seek A/x™ = A/i™(n) — 
A/x'"'(aq) (Fig. [1]), and they express this in terms of the 
free energy change in restricting the coordination state 
of the aqua- ion to n. Thus 

A/i'='^ = A/i'='^(n) - A/i°'^(aq) 

= AmS1+ (n, ANa+ ) - AAi-+ (n, Ak+ ) , (2) 

where Ax is the coordination radius of ion X and n is the 
coordination number under consideration, (n = 8 is used 
for inferring the selectivity of the channel f*|.) Eq. [2]is a 
thermodynamic truism. 



K+(aq) 



A/i'^''(aq) 



Na+(aq) 



A/ig'+(n: Ak+) 



K[H20],|(Ak+ 



A/i<''=(n) 



Na[H20]J(ANa 



FIG. 1: Thermodynamic cycle for calculating A/x''''(n) — 
Aff^{a.q) in terms of A/i5^+ (n, Ak+ ) and A^™ + (n, ANa+ )> 
where A/ix(w,Ax) — /ix('^)Ax) — A'x'(a'<l)- The hydration 
free energy of the ion X in the coordination state n (within a 
coordination sphere of radius Ax) is /iff (n. Ax). The proba- 
bility of observing the n-coordinate state X[H20]n within the 
coordination sphere is 2;x(n, Ax). 

It is in relating A/i™(n, Ax) to the probability of ob- 



2 



serving the n-coordinate state, Xx{n), that the earher 
analysis is deficient. In particular, 



A^°''(aq^ n)' = -fceTln 



X-Na+jn, ANa+) 

xk+ {n, Ak+) 



(3) 



is erroneously equated to A/^'^'^ of Eq. [2] (cf. Eq. 10 in 
Supplementary Information to Ref. 4^1); the prime on 
A/i°''(aq n)' is to note this error. Eq. [3] implies that 
(to within an additivity constant that is the same for 
both K+ and Na+ and without loss of generality taken 
equal to zero), 



A^f=(n,Ax) = -fceTlnxxKAx) 



(4) 



It is rather surprising that the right hand side of Eq. 3] 
is always positive but there is no obvious physical reason 
that the left hand side of Eq. [4]must always be positive. 
Further since J2n 2^x("-, Ax) ~ 1, Eq. [pleads to 



(5) 



If instead of the transition between the aqua-ion to the 
n-coordinate state of the ion, the transition between the 
most probable coordination state of the ion (n) to the 
n-coordinate state is considered, Eq. [3] becomes 



A/.t°'^(n n)' = -feTln 
+ fceTln 



a^Na+ {n, ANa+) 
a:^Na+(^, ANa+) 

XK+(n, Ak+) 
xk+ (^,Ak+)' 



(6) 



where different values of n are implied for Na'^ and K+. 
In the earher study (Fig. 3B, Ref. [1|) the "most popular 
coordinated states" were indicated on the free energy 
landscape (at a free energy value near zero) on which the 
n = 8 coordinate case was also noted. Based on this, 
earlier we ^3,] had concluded that Bostick and Brooks [4] 
had used Eq. [SI when in fact they were using Eq. [31 

Results using Eqs. [3] and [6l Fig. [2 collects the re- 
sults for estimates of the selectivity using Eqs. [3] and [6] 
and the corrected versions of these equations (see be- 
low). Eqs. [3| and [Sj lead to an estimate of selectivity 
about 2 kcal/mol higher than the lower bound of selec- 
tivity (« 3 kcal/mol Q)- Observe that Eq. [6] leads to an 
estimate not much different from that based on Eq. [3] 
This is consistent with the observation Q that in the de- 
scription of the hydration free energy of the ion in terms 
of its various coordination states (and a correction for 
long-range interactions), only coordination states with 
n < h contribute dominantly to the thermodynamics. 
Since for Na"*", n = 6 (ANa+ = 3.1 A), and for K+, n = 7 
(Ak+ = 3.5 A), it is then reasonable to expect that Eqs. [3] 
and [Hi give similar estimates. 



II. DENSITY FLUCTUATIONS AND A^^'=(n,Ax) 

We first note the following identities that readily follow 
from a multi-state generalization of the potential distri- 
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FIG. 2: Estimates of the selectivity with various approaches. 
Red line: the stated lower bound of selectivity Theories 
that lead to an estimate below this line would be regarded 
as inadequate in explaining selectivity. Broken blue and black 
lines: A^i°''(aq 8)' is based on Eq. [3] and Afi'^'^ih 8)' 
is based on Eq. [B] These two estimates differ by about 
0.3 kcal/mol 0.5 ksT). Xj^^+ = 3.1 A and Ak+ = 3.5 A. 
Unbroken blue and black lines: Calculations using Eq. [5]and 
the corrected version of Eq. [6] Symbols: Calculations using 
A = ANa+ — and the correct forms of Eq. [3] (A) and 

Eq. [6](n). Simulation methods are noted in Appendix A. 



bution theorem [1, i, i, 0, S, i, : 

g/3A<5?(aq) ^ ^a;x(n,Ax) •e''^x(".Ax) (7^^) 
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-/3Mx (aq) 



-/3/if=(n,Ax) 



(7b) 



Here p(n. Ax) is the probability of observing n water 
molecules in an ion-free coordination sphere of radius Ax 
in neat liquid water. (We will call the ion-free coordina- 
tion sphere the observation volume.) The set {p(n. Ax)} 
codifies density fluctuations in neat water at the length 
scale of the observation volume [ll|; p(0. Ax) is an essen- 
tial factor in the theory of hydrophobic effects and is a 
measure of the packing interactions involved in hydration 

The following relation was derived earlier 



„/3M5r(aq) ^ ^X{n, Ax) ^g^y(n,Ax) 

pin,Xx) 



(8) 



It can be seen that Eq.[S| satisfies both members of Eqs.[7| 
and also has an intuitively appealing explanation (Fig. [3]). 

Using Eq. [5]in Eq. [21 we obtain 

2^Na+(", ANa+) 



+kBT\n 



a;K+(n, Ak+) 

P(", ANa+) 

p(n,AK+) ■ 



(9) 
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FIG. 3: Schematic description of Eq. [S] The large filled circle is the solute and the small filled circles represent water molecules. 
The dashed circle around the solute denotes the coordination sphere and that in neat water is the observation volume, (a) 
Neat water with an observation volume that is open to the exchange of water molecules with the bulk. The arc with arrow- 
heads on both ends represents the exchangeability of water, (b) p(3) is the probability of observing three water molecules in 
the observation volume in the ensemble of states (a). The free energy change in going from (a) to (b) is —kBTlnp{3). (c) 
The particle is inserted into the observation volume with 3 water molecules. The excess chemical potential of the ion in the 
3-coordinate state is (d) The constraint on the number of water molecules in the coordination sphere is relaxed. The 

probability of state (c) in the ensemble of states (d) is a;(3) and the free energy change in going from (c) to (d) is fcBrina;(3). 
The excess chemical potential of the ion in (d) is /i'^". Thus ^i™ = + fcBrinx(3) - A:BTlnp(3). 



Similarly, an equation for AiJ.'^^{n n) can be obtained. 

Comparing Eq. [3] to Eq. M (or Eqs. H and [5] to Eqs. [8] 
and l7b[ respectively), it is evident that the factor p(n. Ax) 
is omitted in the earlier development This omission 
leaves out physical phenomena of importance in the ther- 
modynamics of hydration, the consequences of which are 
apparent in Fig. [2] (the unprimed A/i™ values). Including 
the role of density fluctuations in neat water at the length 
scale of the observation volume predicts a selectivity less 
than the suggested lower bound of selectivity. Thus con- 
trol of the coordination state of the ion alone is an inad- 
equate explanation of the selectivity of the ion channel. 
As suggested earlier P, 0, Q , acknowledging the chemical 
difference between the coordinating ligand, water in the 
aqueous phase and carbonyl in the S2 site, is necessary 
to understand the selectivity of the ion-channel. 

When A = ANa+ = Ak+ , the second term on the right 
hand side of Eq. [5] is zero. In this instance Eq. [H] is 
the same as Eq. El The estimated A^'='' for A = 3.5 A 
is about 2 kcal/mol and is in fair agreement with the 
1 kcal/mole estimate obtained using Monte Carlo simu- 
lations 3. Most importantly, both these estimates are 
significantly different from the approximately 5 kcal/mol 
obtained using Eq.Owith ANa+ = 3.1 A and Ak+ = 3.5 A 
(Fig. [2) . Thus we can understand the observation that 
when Eq. [3] is used with the same coordination radii, a 
lower selectivity is predicted. The lower selectivity arises 
due to the cancellation of the effects of water occupancy 
in an ion-free coordination volume of the same size for 
both Na+ and K+. 

Before concluding this section, we present some 
thoughts on where the earlier ^] development may have 
faltered. The coordination sphere is an open system 
in contact with an external bath that provides solvent 
molecules at a constant chemical potential, but the coor- 
dination volume and the number of particles it can hold 
are nowhere close to the thermodynamic limit. Thus care 
is required in assigning a thermodynamic potential to 
the partition function of the system. Further, the small- 
system under consideration will be sensitive to fluctua- 
tions. Hence describing the system with different ensem- 



bles may not be equivalent. Equivalence of ensembles is 
strictly valid only in the thermodynamic limit when rela- 
tive fluctuations in the variables conjugate to those used 
to describe the ensemble are small. Thus, perhaps, using 
an open system framework to study a system whose pop- 
ulation is fixed (Cf. Eq. 7 in Supplementary Information 
to Ref. 3) may explain the disagreement with results 
presented here. 

The present development avoids the subtleties noted 
above. The theoretical framework used here [1,[1,0,BB 
[To| is grounded in well-established statistical mechanics, 
and the results (for example. Fig. [3]) agree with physical 
reasoning. 

III. CONCLUSIONS 

The probability of observing a given number of wa- 
ter molecules in neat water in an ion-frec coordination 
volume codifies density fluctuations in the liquid at the 
scale of the coordination volume. Including this effect 
suggests that free energy changes involved in modifying 
the coordination state of an ion in water Q is inadequate 
in explaining the ion selectivity in the ion channel. 

APPENDIX A: METHODS 

NVT simulations of neat water and of Na"*" and K + I13I 
in water were conducted using the NAMD program \W{ . 
Liquid water was described using the TIP3P [lj| model. 
A temperature of 298 K was maintained using a Langevin 
thermostat. The cubic simulation system comprises 306 
water molecules for the pure water simulation; the ion- 
water system comprises an additional ion. The ion is al- 
ways held fixed at the center of the simulation cell. The 
total number density (counting water molecules and the 
ion, if present) is 33.33 nm"'^. For the ion water system, 
after over 10 ns of equilibration data was collected for 
16 ns. Configurations were saved every 100 fs for further 
analysis. For coordination radii from 3.1 A to 3.5 A, 
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the distribution of water molecules in the coordination 
shell of the ion was calculated. An extensively equili- 
brated box of water molecules was further equilibrated 
for 100 ps. Configurations were saved every 250 fs dur- 
ing the course of a production run lasting 1.996 ns. For 
each configuration, observation volumes of radii between 
3.1 A to 3.5 A were centered on 2197 points arranged in 
a cubic lattice within the simulation cell and the distri- 



bution of water molecules calculated. 
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